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Abstract 

Compactness is introduced as a new method to search for the onset of the quark 
matter transition in relativistic heavy ion collisions. That transition supposedly leads 
to stronger compression and higher compactness of the source in coordinate space. 
That effect could be observed via pion interferometry. We propose to measure the 
compactness of the source in the appropriate principal axis frame of the compactness 
tensor in coordinate space. 
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A variety of signatures for the observation of the quark matter transition in relativistic 
heavy ion collisions have been proposed, see e.g. ]I] and references therein. Nowadays it is 
widely believed that irregularities due to such a new state of matter should be seen most 
clearly in excitation functions, i.e. in the bombarding energy dependence, or possibly also in 
the impact parameter or mass-number dependence of various observables. 

However, in spite of circumstantial evidence for such a novel state of matter ||, undis- 
putable "irregularities" in such experimental excitation functions have not been reported 
to date. Most probably a convincing argument that quark matter has indeed been found 
requires several independent measurements of distinct observables pointing to the same con- 
clusion - namely the onset of deconfmement at a particular bombarding energy, impact 
parameter and system size. 

In the present paper we introduce a new observable for the onset of the phase transition. It 
relies on the measurement of the compactness of the source, which is related to the pressure 
and density of the system in the compression and expansion stage of the nucleus-nucleus 
collision. Observables related to the pressure do not only reflect the transient high density 
but rather relate to the order parameter (field), specifying in which thermodynamical state [[| 
the excited matter was created. The compactness can be identified via interferometry: The 
illumination of the baryon source by the pion radiation is subject to experimental scrutinity 
via pion interferometry measurements ||. 

To illustrate the basic idea, let us first discuss compression shocks in heavy-ion collisions 
employing the Rankine Hugoniot Taub shock adiabat (RHTA) solution || |(| of relativistic 
hydrodynamics, 

w>-w > + p(W-^)=o . (1) 

VP Po J 

It relates the energy per baryon in the local rest frame of the compressed matter, W, to 
its baryon density p and pressure P. Wq and po are the energy (=923 MeV) and density 
(=0.16 fm~ 3 ) in the ground state, and the compression factor is p/po- 
Solving for P we obtain 

P = H/ po 7cm ~ 1/7c ; m . (2) 

1 - IcmPo/P 

Simply speaking, the onset of the transition to quark matter at a given incident energy E^g = 
2(7cm 2 _ l)Wo lets the pressure increase less rapidly with W and/or p, and consequently 
higher compression p/po can be achieved than for the case without transition. Now, as 
pV ~ irR\Lp must equal 2A by virtue of baryon number conservation, the longitudinal 
thickness L of the compressed matter is proportional to 1/p. Thus, a transition to quark 
matter leads to a more compact system, just as compact stars with a quark matter core are 
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more compact than pure neutron stars fj§ . Of course, in heavy-ion collisions that expectation 
is based on the behavior of relativistic compression shocks rather than hydrostatic and 
gravitational equilibrium. 

To investigate quantitatively the experimental observable, we perform 3-dimensional (1- 
fluid) calculations of relativistic hydrodynamics for the compression as well as for the sub- 
sequent expansion stage of the reaction. That is, we solve numerically the continuity equa- 
tions for the energy-momentum tensor, d fJi T' lu = 0, and the net baryon current, d^Ng = 0. 
Detailed discussions of (3+l)-d numerical solutions for hydrodynamical compression and 
expansion can be found e.g. in ||. We shall employ two different equations of state (EoS) 
P(W,p): i) a relativistic mean field (RMF) hadron fluid corresponding to baryons and 
antibaryons interacting via exchange of massive scalar and vector bosons, plus free thermal 
pions; the parameters of the lagrangian are fitted to the ground state of infinite nuclear 
matter, in particular the nuclear saturation density, the energy per nucleon, and the incom- 
pressibility. ii) the same EoS as in i) for the low-density phase, supplemented by the Bag 
Model equation of state with bag constant B 1 ^ = 235 MeV for the quark-gluon (QG) phase. 
The phase coexistence region corresponding to this first-order transition is constructed em- 
ploying Gibbs' condition of phase equilibrium, prmf{T, //b) = Pqg(T, //b), where T and \xb 
denote the temperature and the baryon-chemical potential, respectively. For example, for 
p = we find Tq ~ 170 MeV, while at T = phase coexistence sets in at p ~ 4.6po- A more 
detailed discussion of these EoS can be found in 0. 

As an example, we study the compactness in the reactions Au+Au at impact parameter 
b = 3 fm and energy Ef£g = 8A GeV, where the transition to quark matter does occur 
(within the present model). We define the configuration space sphericity tensor, also called 
"compactness tensor" below, as the second moment of the net baryon current. On fixed time 
hypersurfaces we have 

F i:i (t) = J d 3 x Xi Xj N° B (t, x) {pit, x) - Pcut ) . (3) 

We apply an additional density cut p > p cut in the integral to discard spectator matter. 
In the future the cuts and the hypersurface will have to be adapted to the experimental 
conditions. However, that is not crucial for understanding the effect. 

The three eigenvalues /„ are the solutions of the cubic equation det(Fy — fSij) = 0, and 
the eigenvectors e n follow from solving the linear system of equations (Fy — f n bij)e^ n = 0. 
Let e±, e*2 be the eigenvectors defining the reaction plane, e± corresponding to the bigger 
of the eigenvalues f\ and fi- To simplify the discussion we shall assume that the matter 
distribution is symmetric with respect to the reaction plane. 

We can now rotate the coordinate frame around e*3 by an angle cos O = e\ ■ e z , where e z 
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defines the longitudinal (beam-) direction in the lab frame. The rotated compactness tensor 
F* can be written in terms of the eigenvalues f n and orthogonal eigenvectors e* as 



F* = f\e* ® e* + f 2 e; <g> e 2 * + / 3 e 3 * ® e* . 



(4) 



In diagonal form, F* specifies an ellipsoid in configuration space with principal axis along 
e* and radii y/~Jn- For example, cigar patterns, oriented along the z-axis, would correspond 
to fi>f2 = hi e* = e z , e 2 * = e x , e£ = e y . The "compactness" can now be defined as the 
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Figure 1: Contour plot of the time-like component of the net baryon four-current in planes 
transverse to the longitudinal axis e\* of the compactness tensor (Au+Au collision at Ef^ = 
8 AGeV, 6 = 3 fm). Left: Relativistic mean-field EoS without phase transition. Right: with 
phase transition to quark matter. 

Fig. ^] shows typical density distributions in the transverse planes of the coordinate frame 
where F* is diagonal, eq. ([|). The polar angle relative to the lab frame as well as the 
eccentricity 

J3 ~ /2 fr\ 
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are shown as well. As already indicated in the introduction, we find very different eigenvalues 
for the two equations of state, the one without and the other with a first order phase 
transition included. The calculation with transition to quark matter corresponds to higher 
compactness of the baryon distribution. Because the incompressibility dp/ dp or dp/d(Wp) 
is smaller in the presence of the coexistence phase, the compactness tensor is much flatter 
(nearly a factor of two !) in model ii) than in model i). Moreover, our (3+l)-dimensional 
expansion solutions show that the ratio of the in-plane radii J Y2 / fi remains much smaller 
in the case with phase transition, cf. Fig. ||. 
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Figure 2: Ratio of the in-plane radii y/2/ji as it evolves during the compression and ex- 
pansion stages, for the relativistic mean-field EoS without phase transition and for the case 
with transition to quark matter (Au+Au, Ef^ = 8 AGeV, 6 = 3 fm). 

This allows to measure directly the density increase in the high density stage of the 
reaction, if a phase transition occurs. Care must be taken that the investigated range of 
impact parameters constitutes a moderately small bin of centrality values. One should keep 
in mind that the compression factor is affected by the incompressibility evaluated on the 
RHTA (P, but not along a path of fixed specific entropy, because a large amount of entropy 
is being produced in the compression process. 

The predicted change in compactness due to the reduced incompressibility associated with 
a first-order phase transition can be observed experimentally via the novel interferometry 
analysis developed recently by Lisa et at M. The proposed method is quite robust and 



5 



incorporates other interesting information as the configuration space tilt angles []nj, which 
complement the momentum-space flow angles. It avoids cuts in tilted ellipsoids, which are 
not analysed in the appropriate rotated frame, and were the eccentricity and the RMS-radii 
are much less distinct for models i) and ii). 
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